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Abstract 

We present a preliminary measurement of CP- violating asymmetries in fully reconstructed B^^D^*^^Tr^ 
and B^^D^p^ decays in approximately 110 million T(4S') — > BB decays collected with the BABAR 
detector at the PEP-II asymmetric-energy B factory at SLAC. From a maximum likelihood fit 
to the time-dependent decay distributions we obtain for the CP-violating parameters: a^'^ = 
-0.032±0.031 (stat.) ±0.020 (syst.),cj^p = -0.059 ± 0.055 (stat.) ±0.033 (syst.) on the B^^D^tt^ 
sample, a^*'^ = -0.049 ± 0.031 (stat.) ± 0.020 (syst.), cj^*'^ = +0.044 ± 0.054 (stat.) ± 0.033 (syst.) 
on the B°^P>*±7r=F sample, and a^P = -0.005 ± 0.044 (stat.) ± 0.021 (syst.), = -0.147 ± 
0.074 (stat.) ± 0.035 (syst.) on the B^^D^p^ sample. 
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1 INTRODUCTION 



While the measurement of sin 2/3 is now a precision measurement 012]) the constraints on the other 
two angles of the Unitarity Triangle [3], a and 7, are still limited by statistics and/or by theoretical 
uncertainties. This conference paper reports on the measurement of CP-violating asymmetries in 
B'^^D^*^^tt^ and B^^D^ decays jl] in T(4S') — > BE decays, asymmetries which are related 
to I sin(2/3 + 7)1 O El- This analysis updates the results already published in [H] by including a 
new decay mode (B^^D^p^) and a larger data sample (110 instead of 88 million T{4S) — > BB 
decays) . 

The time evolution of B^^D^*^^h'^ decays, where h is a meson made of a and a d quarks, 
is sensitive to 7 because the CKM-favored decay B^^D^*^~^h~ , which amplitude is proportional 
to the CKM matrix elements V^^y*^, and the doubly-CKM-suppressed decay B^^D^*^^h~ , which 
amplitude is proportional to VcdVub interfere due to the B^-B^ mixing. The relative weak phase 
between the two amplitudes is 7, and, when combined with B^B^ mixing, yields a weak phase 
difference of 2/3 + 7 between the interfering amplitudes. 

The decay rate distribution for B^^D^h^ decays is 

f^{rj,At) = X [1 T sm(AmrfAt) ^ rjC cos( Am^At)] , (1) 

where r is the B^ lifetime, neglecting the decay width difference, Am^^ is the B^B^ mixing frequency, 
and At = tree — ^tag is the time of the B^ D^ir^ decay (-Brec) relative to the decay of the other 
B (-Btag)- In this equation the upper (lower) sign refers to the flavor of i?tag as B^ (B^), while 
r] = +1 (—1) and C = + (~) for the final state D~h'^ {D^h~). In the Standard Model, the S and 
C parameters can be expressed as 

2Im(A±) , ^ 1-^^ 

where r = |A+| = 1/|A_| and 

A± = ^A{W^D^t:^)/A{B^^D^t:^) = r±ie-'(2/3+7T5) . (3) 

Here ^ is a function of the elements of the mixing matrix [71, and 6 is the relative strong phase 
between the two contributing amplitudes. In these equations the parameters r and 5 depend on 
the choice of the final state and will be indicated as r^'^ , S^'" , in the B^^D^tt^ case, r^*'", 6^*'^, 
in the B^^D*^7r^ case ^, and r^P, 5^p, in the B^^D^p"^ case. 

Interpreting the S and C parameters in terms of the angles of the Unitarity Triangle requires 
the measurement of the r parameters as detailed in Sec. [21 Since the amplitude in the numerator of 
Eq.|21is suppressed with respect to the one in the denominator, the r parameters are expected to be 
small (~ 0.02) and they cannot therefore be extracted from the measurement of C with the current 
statistics. They can be determined, assuming SU{3) symmetry and neglecting contributions from 
annihilation diagrams, from the ratios of branching fractions B{B^^Di*^^Tr~)/B{B^^D^*'^~Tr~^) 
and B{B^^Di*^^ p')/B{B^^D^*^- p+) ^EM- Since there is no evidence yet of B°^Di*^^p- 
decays, we will report here the result of the measurement of CP-violating parameters defined in 
Eq. but will not be able to include the interpretation in terms of sin(2/3 -|- 7). Nonetheless, the 

^According to Ref. ^ the strong phase is actually 5^ + vr, but this does not affect this measurement. 
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Figure 1: ttt-es distributions in the A£' signal region for, from left to right, the B — > D vr^, 
B'^ D*^iT^ and B^ — > D^p^ sample for events with tagging information. A fit to a Gaussian 
plus a threshold is overlaid. 

addition of the B^ — > D^p^ mode to the analysis is of interest because on one side unexpectedly 
large values for the CP-violating parameters would signal new physics and on the other side when 
the statistics will be sufficient it will allow one more measurement of sin(2/3 + 7), redundant with 
respect to the one made with the other modes. 

2 ANALYSIS OVERVIEW 

This measurement is based on 110 million T(4S') — > BB decays collected with the BABAR detec- 
tor ^Uj at the PEP-II asymmetric-energy B factory at SLAC. We use Monte Carlo simulation of 
the BABAR detector based on GEANT4 [TTj to validate the analysis procedure and to estimate some 
of the backgrounds. The analysis strategy is identical to our previous publication on this topic [S]. 

Candidate B^ D^*^^Tr^ and B^ D^p^ decays are reconstructed with the D*~^ decaying to 
D^Ti^ , where the subsequently decays to one of the four modes K^n^ , K^-n^i:^ , A'^7r+7r^7r+, 
or ivrg-TT"*" 7r~, the decaying into K~it^tt'^ and K^Tr~^, and the decaying into vr'^vr'^. The 
D^*^ candidates are then combined with either a single track or a track and a vr*^ candidate with 
invariant mass in the p window, 620 < mT^-n- < 920MeV/c^. Finally, exploiting the spin properties 
of the decay of a pseudo-scalar meson into another pseudo-scalar and a vector, we require the 
cosine of the angle 9i between the charged pion and the p candidate in the p rest frame to satisfy 
I cos 9h\ > 0.4. 

Signal and background are discriminated by two kinematic variables: the beam-energy sub- 
stituted mass, mES = \/ — P^^, and the difference between the B candidate's measured 
energy and the beam energy, AE = E"^ — (^/s/2). E^ {p%) is the energy (momentum) of the B 
candidate in the e^e~ center-of-mass frame, and ^/s is the total center-of-mass energy. The signal 
region is defined to be \A.E\ < 3a, where the resolution a is mode-dependent and is approximately 
20 MeV as determined from data. Figure Q shows the ttt-es distribution for candidates in the signal 
region. 

To identify the flavor of Stag, each event is assigned by a neural network to one of four hier- 
archical, mutually exclusive tagging categories: one lepton and two kaon categories based on the 
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charges of identified leptons and kaons, and a fourth category for remaining events pp. The effective 
tagging efficiency is (28.1 it 0.7)%. The time difference At is calculated from the measured separa- 
tion along the beam collision axis, Az, between the reconstructed (-Brec) and tagged (-Btag) decay 
vertexes. We determine the B^ec vertex from its charged tracks. The -Btag decay vertex is obtained 
by fitting tracks that do not belong to -Brecj imposing constraints from the B^^c momentum and 
the beam-spot location. The At resolution is approximately 1.1 ps. 

An unbinned likelihood fit is performed to the time distribution of events in this sample. The 
likelihood accounts for possible CP violation on the tag side ^2j: for each tagging category i and 
for each decay mode /i = i^^vr^, D*^7r^ , D"^ 

ft''iv,At) = ^_x[lT(a^T7?&,-r/cf) 

sm{AmdAt) =F V cos{AmdAt)] , (4) 

where in the Standard Model 

a^' = 2r'' sin(2/3 + 7) cos J'' , 
bi = 2r • sin(2/3 + 7) cos (5- , 

= 2 cos (2/3 + 7) (r^ sin 5^ -r^ sin 5-) • (5) 

Here (5-) is, for each tagging category, the effective amplitude (phase) used to parameterize the tag 
side interference. Terms of order r'^ ^ and r-^ have been neglected. Results are quoted only for the 
six and cj^^ parameters, which are independent of the unknowns and (5- (semileptonic B decays 
have no doubly CKM- suppressed contributions and therefore ?'jep=0). The other parameters are 
constrained by the fit, but, as they depend on and 6'^, they do not contribute to the interpretation 
of the result in terms of sin(2/3 + 7). 



3 SAMPLE COMPOSITION 

The background can be separated in two categories, one of which is due to random combinations 
of particles in the event (combinatorial background) and the other is due to B decays into similar 
final states, which therefore has the rriEs distribution similar to the signal (peaking background). 
To separate the combinatorial background, the rriES distribution is fit with the sum of a threshold 
function ^Sj and a Gaussian, with a width of about 2.5MeV/c^, to describe the signal. 



Table 1: Yields, fraction of combinatorial background fcomb and of peaking backgrounds fpeak of 
the selected samples. 



Decay 


yields 


fcomb(%) 


fpeak{%) 


mode 










B^^DV 


7611±97 


8.6 


0.21±0.06 


0.93 ± 0.23 


B^^D*^TT^ 


7068±89 


4.2 


0.13±0.06 


0.93 ± 0.10 


B^^D^p^ 


4400±79 


12.7 


-0.01±0.07 


0.31 ± 0.13 



We considered all possible sources of background peaking in the ttt-es signal region, those coming 
from decays into open-charm final states similar to that of the signal (e.g., B~ D*^7r~,p~ or 
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vr ,p), and those arising from charmless decays into the same final state as the signal. 
We estimated their contributions on MC simulation, varying the branching fractions within errors 
when observed and within the existing upper limits otherwise. The Gaussian yields and the amount 
of peaking background are summarized in Table identified by the source, either neutral or charged 
B mesons. 

In the case of the D*^p^ decays, an additional source of background must be considered, 

which has the same final state, B^ D^tt^tt^, where the tt^tt^ system is not produced by a p 
meson. This background component can be studied by looking at the distribution of m^7r and 
cos Oh- When the vr^ and vr" come from a p meson, rnvrTr follows the p lineshape, while cosO^ is 
distributed as {cosOh)'^- In order to satisfy the Bose-Einstein statistics, vr^ and vr'' must come from 
a resonance with odd spins (J = 1,3,...), and therefore, to first approximation, we only consider 
excited states of the p meson. In the mass range of interest, p(1450) is the only possible candidate. 
For B^ — > D^p(1450) decays, the m,r7r distribution would be peaked at higher masses (the />(1450) 
pole mass is (1465 ± 25) MeV/c^, its width is (400 ±60) jTi), but would have the same (cos^/j)^ 
distribution. In order to have a cos 9h distribution different from the signal, the background would 
have to come from non-resonant B^ D^n^ir^ decays. We therefore consider the possibility of 
having three components: the signal, B^ D^p^, B^ — > Z?^p(1450), and an S-wave non-resonant 

component, 5° D'^{TT'^TT^)nr- 

The presence of these two additional components would imply that the value of 
A{B^^D^n^'ir°)/A{B^^D^TT^iT^) (and thus A^'', see Eq.^ is not constant with m^^: 

A(B°^L>^7r±7r°)/A(5°^Z)^7r±7r°) = (6) 
A{BO ^ D- p+)Ap{m^^) + A{BO ^ Z)-p(1450)+)A^(i450)("i--) + ^(^° ^ D' 
A{BO ^ D-p+)Ap{m^^) + A{BO ^ Z)-p(1450)+)^p(i450) ("i--) + ^ D- 

Here Ap(m7r7r)5 ^p(i450) ("^tttt), and A„r("T.7r7r) are the amplitudes of the three components and de- 
pend on rriTrn'^ . 

FigurelSlshows the comparison of the m-j^Tr and cos 9h distributions between data and a simulation 
of pure B^^D^p"^ decays. The good agreement suggests that the contribution of higher resonances 
is negligible. To quantify this statement, we perform a fit to the mj^T^ distribution including the 
lineshapes of the three components, and extract the relative amplitudes and phases. Given this 
set of amplitudes and phases, we generate the m,r7r distribution and, for each value of m^^, the 
appropriate At distribution. We then fit these samples with the likelihood described in Eq.@, 
which ignores the dependence of on m-^T^. We repeat this procedure for several sets of amplitudes 
and phases according to the measured covariance matrix and find the bias induced to the a (c) 
parameters to be at most 0.0018 (0.0047). This maximum bias will be included as a systematic 
error. 

An unbinned maximum likelihood fit is performed on the selected B candidates using the signal 
At distribution in Eq.(|3I), convoluted with a three-Gaussian resolution function. The probabil- 
ities of incorrect tagging (wi) are accounted for by multiplying the a^, cf^ parameters and the 
cos(AmdAt) term by the dilutions Di = 1 — 2wi. The resolution function and the tagging param- 
eters are consistent within errors with previous BABAR analyses 

The combinatorial background is parametrized as the sum of a component with zero lifetime and 
one with an effective lifetime fixed to the value obtained from simulation. The fraction of each back- 
ground component is determined from the events in the uies sidebands, 5.2 < tties < 5.27GeV/c^, 

^In a Dalitz analysis they would depend also on the D vr invariant mass, but here we integrate over it to study a 
single variable at a time. 
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Figure 2: Sideband-subtracted 771^^ and cos Oh distributions on data (dots) and B^^D^p^ simu- 
lation (open histogram). 



while the At resolution is a double- Gaussian fitted on the same data as the CP-violating parame- 
ters. The charmed background coming from mesons is modeled by an exponential with the B^ 
lifetime, and its amount is fixed to the value predicted by simulation. The charmed and charmless 
backgrounds from B^ mesons are neglected in the nominal fit, but are considered in evaluating the 
systematic uncertainties. 



4 RESULTS 

From the unbinned maximum likelihood fit we obtain: 

^Dtt ^ -0.032 ± 0.031 (stat.) , cg^ = -0.059 ± 0.055 (stat.) 

^D-TT ^ -0.049 ±0.031 (stat.) , cg'^ = 0.044 ± 0.054 (stat.) 

a^P = -0.005 ±0.044 (stat.) , c& = -0.147 ± 0.074 (stat.) 



Table m shows the contributions to the systematic uncertainty. Errors are found to be indepen- 
dent of the Brec reconstruction mode because the systematic effects are dominated by uncertainties 
in the Btag reconstruction. In the table we compare the results obtained with B'^^D^*^^tt'^ decays 
and B^^D^p^. 

The impact of a possible mismeasurement of At {(tm) has been estimated by comparing the 
different parameterizations of the resolution function, varying the position of the beam spot, and 
the absolute z scale within their uncertainties, and loosening and tightening the quality criteria on 
the reconstructed vertex. We also estimate the impact of the uncertainties on the alignment of the 
silicon detector (SVT) by repeating the measurement on simulated events, intentionally misaligning 
the SVT in the simulation. As systematic uncertainty of the fit (fifit), we quote the upper limit 
on the bias on the and c^, as estimated from samples of fully simulated events. The model 



12 



Table 2: Systematic uncertainties on the a and c parameters. 









B^^D^p^ 


Source 










Vertexing {a At) 


0.015 


0.026 


0.017 


0.031 


Fit (cjfit) 


0.011 


0.019 


0.009 


0.016 


Model ((Tmod) 


0.006 


0.007 


0.0007 


0.0015 


Tagging (crtag) 


0.004 


0.0034 


0.0028 


0.0033 


Background (o"bkg) 


0.0012 


0.0027 


0.006 


0.0031 


Dependence from m^^o (cAdep) 






0.0018 


0.0047 


Total (cJtot) 


0.020 


0.033 


0.021 


0.035 



error (fTmod) contains the uncertainty on the B^ lifetime and Am^, varied by the uncertainties on 
the world averages 0, and the impact of neglecting higher order terms in r or r- in Eq.@. The 
tagging error ((Jtag) is estimated considering possible differences in tagging efficiency between B^ 
and B^ and allowing for different At resolutions for correctly and incorrectly tagged events. We also 
account for uncertainties on the background (cTbkg) by varying the effective lifetimes, dilutions, tues 
shape parameters, signal fractions, and background CP asymmetry up to five times the expected 
CP asymmetry for signal. For the B Dp decay we also include the maximum bias of the a and 
ciep parameters due to the possible dependence of A on the vrvr^ invariant mass (fJAdep)) as discussed 
in sectional 

5 CONCLUSIONS 

We studied the time evolution of fully reconstructed B^^D^*^^Tr^ and B^^D'^ p^ decays in a 
data sample of 110 million T(45') — > BB decays. CP-violation arising from the interference of 
the CKM-suppressed and the CKM-favored amplitudes is expected to be small but sensitive to 
sin(2/3 + 7). 

The CP-violating parameters defined in Eq. Elare measured to be 



^Dtt ^ -0.032 ± 0.031 (stat.) ± 0.020 (syst.) 
^D*n ^ -0.049 ± 0.031 (stat.) ± 0.020 (syst.) 
a^P = -0.005 ±0.044 (stat.) ±0.021 (syst.) 

No significant CP asymmetry is observed thus far. 
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